A&A, 686, A280 (2024) 
https://doi.org/10.1051/0004-6361/202348242 
© The Authors 2024 


tronomy 
Astrophysics 


A too-many dwarf satellite galaxies problem in the MATLAS 
low-to-moderate density fields 


Kosuke Jamie Kanehisa!:2®, Marcel S. Pawlowski!®, Nick Heesters*®, and Oliver Müller’ © 


' Leibniz-Institut für Astrophysik Potsdam (AIP), An der Sternwarte 16, 14482 Potsdam, Germany 
e-mail: kkanehisa@aip.de 
2 Institut fiir Physik und Astronomie, Universitat Potsdam, Karl-Liebknecht-StraBe 24/25, 14476 Potsdam, Germany 
3 Institute of Physics, Laboratory of Astrophysics, Ecole Polytechnique Fédérale de Lausanne (EPFL), 1290 Sauverny, Switzerland 


Received 11 October 2023 / Accepted 11 April 2024 


ABSTRACT 


Context. Dwarf galaxy abundances can serve as discernment tests for models of structure formation. Previous small-scale tensions be- 
tween observations and dark matter-only cosmological simulations may have been resolved with the inclusion of baryonic processes; 
however, these successes have been largely concentrated on the Local Group dwarfs the feedback models were initially calibrated on. 
Aims. We investigate whether the ACDM model can reliably reproduce dwarf abundances in the MATLAS low-to-moderate density 
fields that are centred upon early-type host galaxies beyond the Local Volume. 

Methods. We carried out mock observations of MATLAS-like fields with the high-resolution hydrodynamic simulation IlustrisTNG- 
50. We used matching selection criteria and compared the properties of dwarfs contained within them with their MATLAS analogues. 
Results. Although simulated MATLAS-like dwarfs demonstrate photometric properties that are consistent with the observed galaxy 
population and follow the same scaling relations, TNG50 underestimates the number of dwarf galaxies in isolated MATLAS fields at 
the 60 level. This significance maintained within crowded fields containing more than a single bright host. Our 55-62% estimate of 
the fraction of background galaxies is in agreement with estimates by MATLAS, but is wholly insufficient to alleviate this discrepancy 
in dwarf abundances. Any incompleteness in the observed fields further exacerbates this tension. 

Conclusions. We identified a “too-many-satellites” problem in ACDM, emphasising the need for the continued testing and refining 


of current models of galaxy formation in environments beyond the Local Group. 
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1. Introduction 


The concordance ACDM model has been largely success- 
ful in matching astrophysical observations at cosmological 
scales, such as the accelerated expansion of the Universe 
(Planck Collaboration XIII 2016), the CMB power spectrum 
(Spergel et al. 2007), and filamentary large-scale structure 
(Tempel et al. 2014). However, comparing the expectations from 
cosmological simulations with observed galaxies and their dwarf 
satellites has revealed a number of tensions at galactic scales 
(Bullock & Boylan-Kolchin 2017; Sales et al. 2022). 
Abundances of dwarf galaxies are particularly sensitive to 
the feedback prescription adopted. In the ‘missing satellites’ 
problem, dark matter-only simulations predict many more sub- 
halos around Milky Way and M 31-mass galaxies than what 
has been observed with respect to dwarf satellites (Klypin et al. 
1999; Moore et al. 1999). This discrepancy arises from the 
quenching of galaxy formation in low-mass subhalos from 
gas heating due to the UV background from cosmic reionisa- 
tion (Bovill & Ricotti 2009; Revaz & Jablonka 2018). Modern 
hydrodynamical and high-resolution ACDM simulations consis- 
tently match the luminosity function of the Milky Way and M 31 
dwarfs; thus, there no longer appears to be a “missing satellites’ 
problem within the Local Group (Engler et al. 2021). 
Conversely, we are faced with a ‘too-big-to-fail’ (TBTF) 
problem (Boylan-Kolchin et al. 2011) with respect to the high- 
mass end of the subhalo mass function. The dark matter masses 


of the brightest Milky Way satellites, as inferred from their inter- 
nal velocities, are significantly lower than the most massive sub- 
halos around simulated analogues (Boylan-Kolchin et al. 2012). 
These dark subhalos are sufficiently massive to enable the effec- 
tive cooling of their gas content and, thus, they are too big to fail 
in forming visible galaxies. Yet, with baryons, the increased host 
galaxy potential enhances the tidal stripping and mass loss of 
their satellites (Brooks & Zolotov 2014), while the formation of 
cored profiles effectively reduce the host halo mass and, accord- 
ingly, the expected subhalo abundance (Peñarrubia et al. 2010). 
There is now a consensus among modern hydrodynamic simula- 
tions that there is no TBTF problem for Milky Way and M 31- 
mass hosts (Sawala et al. 2016; Samuel et al. 2020); however, 
we refer to Pawlowski et al. (2015) for an alternative viewpoint. 

Since the baryonic solutions implemented in cosmologi- 
cal simulations are calibrated for Local Group-like environ- 
ments, we highlight the need to verify whether these successes 
remain robust beyond our cosmic neighbourhood. For instance, 
Smercina et al. (2018) searched for dwarfs within 150 kpc of 
the Milky Way-mass spiral galaxy M 94 (D = 4.2Mpc; 
Radburn-Smith et al. 2011) and found only two ‘classical’ satel- 
lites with M, > 4 x 10° Mo. The four other comparably massive 
hosts with complete ‘classical’ samples host an average of 6—12 
such satellites, and only <0.2% of M 94 analogs in the EAGLE 
simulation (Schaye et al. 2015) host similarly sparse satellite 
systems. Similarly, a study of seven low-mass spiral galaxies in 
the galactic neighbourhood (Miiller & Jerjen 2020) found that 
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high-resolution dark matter simulations appear to systematically 
overestimate the observed satellite abundances (although this 
discrepancy may arise from observational biases). Conversely, 
the satellite luminosity function around Centaurus A — a massive 
elliptical galaxy at 3.8 Mpc that has been the target of multiple 
deep surveys — is shown to be consistent with expectations from 
the hydrodynamic IlustrisTNG-100 simulation to within a 90% 
confidence interval (Miiller et al. 2019). 

Recently, Miiller et al. (2024) studied the abundance of 
dwarf galaxies around M 83, a barred spiral at D = 4.9Mpc 
(Jacobs et al. 2009). with a mass (and, thus, the expected satel- 
lite abundances) that is comparable to the Milky Way values. 
M 83 has 13 confirmed satellites within a projected separation of 
330 kpc down to a limiting magnitude of My = —10 mag. M 83 
analogs in the high-resolution HlustrisTNG-50 simulation, how- 
ever, demonstrate sparse satellite distributions with luminosity 
functions inconsistent with that of the M 83 system at a 30 confi- 
dence level — thus constituting the so-called ‘too-many-satellites’ 
problem. The actual degree of tension is likely underestimated 
due to incompleteness in the observed satellite sample. 

These ambiguous results suggest that models of galaxy for- 
mation in ACDM may still struggle to consistently predict the 
abundances of dwarf satellites around host galaxies beyond the 
Local Group. And yet, the works above only consider dwarf 
abundances of a few individual systems, which may conceiv- 
ably be over- or under-populated as a result of some unique 
evolutionary history. In this paper, we study dwarf galaxy abun- 
dances within the low-to-moderate density fields in the MAT- 
LAS survey (Duc et al. 2015; Habas et al. 2020; Bilek et al. 
2020; Poulain et al. 2021), across a total of 150 fields centred 
upon early-type host galaxies beyond the Local Volume. We 
then mock-observe analogous fields using identical photomet- 
ric selection criteria in the state-of-the-art hydrodynamic simu- 
lation IllustrisTNG-50, which features a resolution rivalling or 
even exceeding zoom-in simulations which consistently predict 
dwarf properties to the lower end of the subhalo mass function. 
In this way, we tested for the first time whether dwarf abun- 
dances in MATLAS are consistent with expectations from con- 
cordance cosmology. 


2. Method 
2.1. ATLAS? and the MATLAS survey 


The Mass Assembly of early-Type gaLAxies with their fine 
Structures (MATLAS) deep imaging survey targets a volume- 
complete sample of early-type galaxies (ETGs) acquired from 
the ATLAS? legacy program (Cappellari et al. 2011). ATLAS?? 
consists of 260 elliptical and lenticular galaxies with K-band abso- 
lute magnitudes of Mx < —21.5, selected from a parent sample 
of 871 massive galaxies. The targeted ETGs are located at decli- 
nations of |6 — 29°| < 35° and Galactic latitudes of b > 15°, and 
lie within 10—45 Mpc (z < 0.01) of the Milky Way. Together with 
the Next Generation Virgo cluster Survey (NGVS; Ferrarese et al. 
2012), MATLAS conducted deep optical imaging of the regions 
around the ATLAS?” ETGs. NGVS mapped fields around 58 of 
the ATLAS? targets between 2009 and 2013, while MATLAS 
has since imaged most of the remaining ETGs. 

The MATLAS survey consists of a total of 150 fields around 
180 ETGs and 59 late-type galaxies (LTGs) from ATLAS*”’s 
parent catalogue, each with dimensions of 63’ x 69’ and a res- 
olution of 0.187 arcsec pix~'. Multiband observations were per- 
formed using the MegaCam on the 3.6-m Canada-France-Hawaii 
Telescope (CFHT) between 2010 and 2015 (Duc et al. 2015). 
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MATLAS was designed to study low surface brightness features 
in the outskirts of the ATLAS*> ETGs to a local depth in the 
g-band of u4 = 28.5—-29 mag arcsec™? (Duc et al. 2014); there- 
fore, it is ideal for resolving populations of faint and low-surface 
brightness dwarf galaxies around their massive hosts. 


2.2. MATLAS dwarf catalogue 


MATLAS adopts both a fully visual and semi-automatic 
approach to generate its dwarf catalogue. The visual approach is 
used purely to calibrate the dwarf selection criteria used by the 
latter and the authors identified 1349 dwarfs across all 150 fields. 
The semi-automatic catalogue first uses SOURCE EXTRACTOR 
(hereafter SEXTRACTOR; Bertin & Arnouts 1996) to generate an 
automated list of dwarf candidates, followed up by a multi- 
stage visual cleaning procedure. It was found that a linear cut 
in average surface brightness, (u,), and apparent magnitude, 
mg, resulted in the cleanest separation between the visually 
selected dwarfs and non-dwarf detections from a representative 
field around NGC 4382 (Habas et al. 2020). Substructures could 
only be identified in more massive galaxies — hence confirm- 
ing their background nature — when their area was larger than 
75 pix? (23.5 arcsec’) and all smaller objects were removed from 
the selection. The catalogue was then cleaned manually by the 
authors and dwarfs were classified by their visual morphology. 

The final dwarf catalogue contains 2210 unique dwarf can- 
didates and is dominated by dwarf ellipticals (dE, 73.4%), while 
the remaining 26.6% are classified as dwarf irregulars (dIrr). One 
of the 150 MATLAS fields is devoid of detectable dwarfs due 
to intense light contamination from a nearby star. The remain- 
ing 149 fields host between 2 and 79 dwarf candidates per field. 
Furthermore, 67 fields contain no bright galaxies (defined by a 
magnitude of <Mx + 1), other than the field’s target ETG, while 
a further 22 fields lack bright galaxies in close projection to the 
central ETG within a square 175 kpc region. 

Dwarf galaxies are generally assumed to be satellites of the 
central ETGs in their respective fields. There are 325 dwarf 
candidates (15%) with at least one independent distance esti- 
mate, spanning a range between 5.1 and 100.6 Mpc. Further- 
more, 90% of these dwarfs lie within the 10—45 Mpc region 
populated by the ATLAS?” ETGs. Intrinsic properties includ- 
ing effective radius, Re, Sersic index, n, and axis ratio, b/a, were 
derived using GALFIT (Peng et al. 2010) in Poulain et al. (2021). 
Since a majority of the MATLAS dwarfs lack associated dis- 
tances, their properties were calculated under the assumption 
that each dwarf is located at the same distance as their field’s 
targeted host galaxy. The authors confirm the dwarf nature of 
the catalogued objects using a cut of M, > —18, and 80% of 
dwarfs appear to be physically associated to their field’s cen- 
tral ETG with a relative radial velocity within 500 km s™!. Of 56 
dwarfs with additional MUSE spectroscopy, 42 (75%) demon- 
strate radial velocities within 1000-3000 km s™!, thus consistent 
with the ATLAS? galaxies (Heesters et al. 2023). 


2.3. Host galaxy selection in TNG50 


We use data from the IllustrisTNG suite! of large-volume hydro- 
dynamic cosmological simulations (Pillepich et al. 2018). All 
runs contained therein adopt cosmological parameters taken 
from Planck (Planck Collaboration XII 2016): Qa = 0.6911, 
Om = 0.3089, and h = 0.6774. We specifically used the highest- 
resolution TNGSO run (Lbox = 51.7 Mpc) which (with a bary- 
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Fig. 1. Properties of ATLAS? host galaxies and their TNGSO analogues. Left: luminosity functions of TNG50 and ATLAS?” hosts are plotted in 
solid and dotted lines, respectively. The latter corresponds to Fig. 3 in Cappellari et al. (2011). Samples containing only early-type and late-type 
centrals are denoted in red and blue respectively. Middle: mass distribution of the dark matter halos hosting the TNGSO analogues. Right: distance 
distribution of ATLAS* targets, smoothened via a Gaussian kernel for re-sampling in this work. The MATLAS survey targets a subset of the 


early-type ATLAS? hosts (solid red line). 


onic particle mass of mga, = 8.5 x 10* Mo) best resolves popu- 
lations of MATLAS-like dwarfs with estimated stellar masses of 
5.8 < log;9(M./Mo) < 9 (see Habas et al. 2020). 

IllustrisTNG adopts the classic friends-of-friends (FOF) 
approach to identify virialised dark halos, then populating them 
with subhalos using SUBFIND (Springel et al. 2001). To select 
ATLAS? -like host galaxies in TNG50, we first identified all 
subhalos considered by SUBFIND to be the central subhalo of 
some parent halo. ATLAS??-like hosts are selected by mirroring 
Cappellari et al. (2011)’s criterion and adopting a K-band abso- 
lute magnitude cut of Mg < —21.5. Since TNGSO contains sev- 
eral highly luminous galaxies up to Mx ~ —28 that lack analogues 
in the parent ATLAS?» catalogue, we imposed an additional con- 
straint of Mx > —26.0 and recover a total of 925 TNG hosts. 

We classified host galaxies by their Hubble morphology into 
early-types (ETGs) and late-types (LTGs). This distinction is 
made by referring to Zana et al. (2022)’s kinematic decomposi- 
tion of M, > 10° Mo subhalos in TNGSO into five distinct com- 
ponents — the thin disc, thick disc, pseudo-bulge, bulge, and stellar 
halo. While the original work defines a given galaxy’s disc-to-total 
mass ratio D/T as the combined mass fraction of the former three 
rotating components, we argue that only the thin and thick disc 
would be observationally distinguishable; thus, we excluded the 
pseudo-bulge from our D/T calculation. Since there is no clear 
distinction between ETGs and LTGs in this metric, we took advan- 
tage of ATLAS?’s volume-complete nature and assume the frac- 
tion of observed ETGs (260/871 = 0.299) also holds in TNGSO. 
The 30th percentile of the full D/T distribution at 0.575 is hence 
set as the threshold between early-types and late-types, and we 
identify 287 ATLAS*»-like ETGs. 

The properties of the sampled ATLAS??-like hosts in 
TNGSO are shown in Fig. 1. The simulated luminosity function 
roughly agrees with that in ATLAS?” (Cappellari et al. 2011), 
although TNG50 demonstrates an excess of luminous early-type 
hosts at Mx < —24. The sample of TNG50 analogues appears to 
be dominated by Milky-Way and M31-mass objects within halos 
of 10!!5 Mo < Ma0 < 10!25 Mo. 


2.4. Mock-observing MATLAS-like fields 


We mock-observed 63’ x 69’ fields centred upon the selected 
TNG hosts to resemble the MATLAS low-to-moderate den- 


sity fields as closely as possible. ATLAS*’s distribution of the 
180 ETGs targeted by MATLAS (Cappellari et al. 2011) were 
smoothed using a Gaussian kernel. The obtained kernel density 
estimate (KDE) was truncated at 10 Mpc and 45 Mpc to mirror 
the observed sample (see the right-hand panel in Fig. 1). We 
generate Nea) = 10 realisations for each early-type TNG host 
galaxy, each with a mock-observation distance, Dhost, indepen- 
dently drawn from the ATLAS? KDE. 

For each realisation, we drew an observing direction, 7, from 
an isotropic distribution and define an associated cone with an 
opening angle of Ofeia = 46.7 arcmin (thus fully containing the 
final 63’ x 69’ survey footprint centred upon the target host) and 
observation depth, Dops, aligned with 7. An observer is placed at 
a distance of Dos from each ATLAS*”-like host galaxy along 
—? and we proceeded to mock-observe all luminous subhalos 
within a 63’ x 69’ field around the host across multiple periodic 
volumes until D,ps = 100 Mpc was reached. 

Since it is unfeasible to  straightforwardly replicate 
Habas et al. (2020)’s visual inspection of dwarf candidates for 
the mock-observed subhalos, we focussed on adopting the quan- 
titative selection criteria used to compile the final MATLAS 
dwarf sample. The greatest challenge lies in estimating effec- 
tive (half-light) radii Re for TNG’s luminous subhalos, a rigor- 
ous approach to which would require the individual modelling 
of Sérsic profiles for each subhalo while accounting for their rel- 
ative orientation and morphology. We instead approximated R, 
with a straightforward conversion from stellar half-mass radii 
Rı/2, which is provided by TNG’s Subfind group catalogues. 
Price et al. (2022) found that the ratio between R1/2 and Re is 
generally consistent over varying Sérsic indices and, rather, it 
primarily depends on the galaxy’s projected axis ratio, go. The 
full population of MATLAS dwarfs (Habas et al. 2020) demon- 
strates a distribution of go between 0.3 and 1.0, with a median 
value of roughly qo = 0.7. The range of MATLAS qo from 
GALFIT covers 1.05 < Rjj2/R. < 1.3, and we adopt a cor- 
responding static R1/2/Re ratio of 1.2 (see Fig. 2 in Price et al. 
2022). 

The apparent g-band magnitude of TNG’s luminous subhalos 
is derived from their absolute magnitude M, in the usual manner 
using: 


mg = M; + Slog, 9(D/10 pe), (1) 
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when mock-observed at a distance D, although we do not 
account for the effects of dust attenuation. Each subhalo’s aver- 
age surface brightness (u4) within Re (in units of mag arcsec”) 
is calculated from mg and the angular expression of Re as: 


(Ug) = 0.75 + mg + 2.5logio(x (R./arcsec)”), (2) 


where the factor 2.5 log;g(2) = 0.75 arises from considering the 
flux observed only within Re. 


2.5. Selection criteria for MATLAS-like dwarfs 


The distribution of observed MATLAS dwarfs in (u,)—m, space 
is plotted in the upper left panel of Fig. 2. We differentiate 
between dwarfs and more massive galaxies by implementing 
Habas et al. (2020)’s linear photometric cut of 


(Lg) = 0.32m; + 18.5. (3) 


While the authors calibrate this cut for each field to best separate 
the low-surface-brightness region of the main locus of galaxies 
with the visually confirmed MATLAS dwarfs, these shifts gen- 
erally remain within (u4) offsets of less than 0.1 mag arcsec”. 
Hence, we adopted Eq. (3) (blue dashed line) as a static 
threshold for the remainder of this work. While MATLAS is 
capable of detecting substructures to local depths of uw, = 
28.5—29 mag arcsec™*, the lowest average surface brightness 
within Re recorded within the final dwarf catalogue is only 
26.7 mag arcsec”. Hence, we defined a criterion for average sur- 
face brightness of: 


(Hg) < 26.7 mag arcsec *, (4) 


for dwarfs in the mock-observed TNG fields (green dashed line). 

The region in (u,)—m, space encompassed by the two cuts 
defined above contains a population of visually bright galaxies 
at m, < 15 that is not present in the MATLAS catalogue. Their 
absence is likely due to the visual inspection by Habas et al. 
(2020) correctly identifying them as non-dwarfs. Since all MAT- 
LAS dwarfs demonstrate GALFIT and SEXTRACTOR magnitudes 
of mg > 14.86, we corrected for this discrepancy by imposing an 
additional, conservative cut at: 


my > 15, (5) 


(red dotted line). Finally, MATLAS rejects all dwarf candidates 
with an ISOAREA parameter in SEXTRACTOR of 75 pix”. While 
ISOAREA has no exact physical analogue available for TNG sub- 
halos, we find that this rejection cut corresponds to a minimum 
angular effective radius of: 


Re > 2 arcsec, (6) 


(black dotted line) in the MATLAS catalogue, which we adopted 
as our last selection criteria for MATLAS-like dwarfs. This 
cut in R, serves to further remove background interlopers from 
our sample, as only <0.1% of physically bound MATLAS-like 
dwarfs are rejected in this final step. 

The population of galaxies in our mock-observed TNG fields 
consist of satellite galaxies of the ETGs targeted by each field, 
as well as a fraction of interlopers — foreground or background 
objects that are not physically associated with their field’s cen- 
tral galaxy. An accurate distinction can be made by whether 
a MATLAS-like dwarf is considered to be bound to their pre- 
sumed host galaxy’s halo in SUBFIND’s halo catalogue. A more 
observationally analogous approach uses the difference in radial 
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velocities between the given dwarf and its field’s central ETG 
while taking their recession velocities due to the Hubble flow 
into account. Habas et al. (2020) estimated that using a cut of 
|Av| < 500km s™!, around 80% of dwarf candidates with inde- 
pendent distance estimates appear to be satellites of the ETG or 
LTG with the smallest angular separation. We adopted the same 
threshold below to test whether we obtain a comparable fraction 
of physically associated dwarfs in TNG’s MATLAS-like fields. 


3. Results and discussion 
3.1. Properties of MATLAS-like dwarfs 


Mock-observed properties of TNG50 subhalos are plotted in 
the left-hand panels of Fig. 2. The distribution of all luminous 
TNG subhalos bound to their field’s central host galaxy follows 
a roughly linear relation in (44,)—m,. Dwarfs considered to be 
MATLAS-like lie in the region enclosed by the four selection 
criteria in Eqs. (3)—(6) and exhibit stellar masses of the order 
of 107—108 Mo. The observed MATLAS dwarfs in the upper 
panel inhabit a region consistent with our mock-observed ana- 
logues, but extend to higher mg for a given (ug). This shift can 
be accounted for if a fraction of the MATLAS dwarfs are back- 
ground interlopers unassociated with their supposed host galax- 
ies. In Fig. 2, we demonstrate that these outlying MATLAS 
dwarfs are consistent with a background distribution at distances 
up to 90-100 Mpc. This maximum is in excellent agreement 
with the 13% of MATLAS dwarfs with independent distance 
estimates — which cover a range of distances up to 100.6 Mpc — 
and justifies our choice of observation depth of Dobs = 100 Mpc 
in Sect. 2.4. 

The right-hand panels of Fig. 2 shows intrinsic properties 
of all luminous TNG subhalos. MATLAS-like dwarfs (plotted 
as white circles) appear to follow the standard scaling relation 
between M} and Re. The full distribution of observed MAT- 
LAS dwarfs partially lies along the region inhabited by TNG’s 
MATLAS-like analogues, but also appreciably extends to larger 
M, and lower Re. Since 85% of MATLAS dwarfs lack distance 
estimates, both Re and M} were calculated by assuming they 
are located at the same distance as their supposed host galaxy 
(Habas et al. 2020). When making this assumption for the lumi- 
nous TNG subhalos, they span a much larger region in M,—R, 
space that is now consistent with the MATLAS dwarfs. The 
fact that the MATLAS dwarfs also extend beyond the (ug) = 
26.7 mag arcsec~” threshold (despite following this cut in the 
left-hand panels) also hints at the unreliability of their size and 
absolute magnitude estimates. The few MATLAS dwarfs with 
distance estimates (blue triangles) demonstrate a much closer 
alignment with TNG’s MATLAS-like dwarfs and generally do 
not extend beyond the (ug) threshold. Overall, the MATLAS- 
like dwarfs mock-observed in TNG50 display photometric and 
intrinsic properties that are fully consistent with dwarfs observed 
in the MATLAS fields, especially when assuming a degree of 
contamination by background interlopers. 


3.2. MATLAS-like dwarf populations 


The number of dwarf galaxies contained within the MATLAS 
fields depends heavily on the central host galaxy’s degree of iso- 
lation. Crowded fields containing multiple ETGs and/or LTGs 
have enhanced dwarf counts due to the stacking of multiple dis- 
tinct satellite populations. To control for this bias, we focus on 
fully isolated fields wherein the central ETG has no other bright 
galaxy with a magnitude of <Mx + 1 within a 3D distance of 
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Fig. 2. Scaling relations for MATLAS dwarfs and luminous subhaloes in the TNGSO simulation. Left: selection criteria for MATLAS-like dwarfs 
mock-observed in TNG50. The blue line indicates the linear cut between regular galaxies and dwarfs adopted in MATLAS (Habas et al. 2020). 
The green, red, and black lines respectively represent cuts on average surface brightness (u4), apparent magnitude m,, and apparent effective radius 
motivated by the distribution of the MATLAS dwarfs (see upper panel). The observed dwarfs demonstrate photometric properties consistent with 
their TNG counterparts, except for a population at larger m, than expected in the simulated (u,)—m, relation. This discrepancy can be explained 
by assuming that a fraction of the MATLAS dwarfs are background galaxies with distances up to 100 Mpc. Right: physical properties of luminous 
subhalos in TNG’s mock-observed fields. MATLAS-like dwarfs that satisfy the selection criteria are indicated as white circles. Distances are not 
available for around 85% of the MATLAS dwarfs, and they are assumed to lie at the same distance as their field’s target host galaxy. The TNG 
dwarfs, when making this assumption, also follow a similar extended distribution in Re— M, space. MATLAS dwarfs with known distances follow 


the true TNG distribution more closely. 


1 Mpc — a threshold selected to distinguish isolated hosts from 
the Local Group pair ~800 kpc apart, see Appendix A. We disre- 
garded galaxies without apparent k-band magnitudes since they 
are likely to lie on the fainter, low-mass end of the luminosity 
function. We searched for potential companion galaxies in the 
Heraklion Extragalactic Catalogue (HECATE; Kovlakas et al. 
2021), a value-added all-sky galaxy catalogue containing over 
200000 galaxies within 200 Mpc. Thus, we identified 48 fully 
isolated MATLAS fields that satisfy these criteria. 

The left-hand panel in Fig. 3 plots the cumulative distribu- 
tion of dwarf populations in the isolated MATLAS fields (solid 
red line). We identified isolated MATLAS-like fields in TNG50 
using the same criteria and plot their dwarf counts in black. In 
stark contrast to the MATLAS median of ten dwarfs per isolated 
field, the mock-observed TNG fields only contain a median pop- 
ulation of three MATLAS-like dwarf galaxies. A two-sample 
Kolmogorov—Smirnov (KS) test reveals that the two distribu- 
tions are inconsistent with being drawn from the same parent dis- 
tribution at a 6.70 confidence level (D = 0.503, p = 2.9x 1071). 
This discrepancy only widens when solely considering TNG 


dwarfs with radial velocities within 500 kms~' of their field’s 
central ETG; thus, MATLAS-like fields in TNG50 are found 
to only contain a median of one such dwarf. When restrict- 
ing our search to dwarfs bound to their field’s host galaxy, the 
median satellite population drops to zero. While we focus on iso- 
lated fields for this analysis, it is evident (see Fig. 3) that more 
crowded fields in TNG and MATLAS demonstrate a similar or 
greater degree of tension. We also note that this result is robust 
to our adopted threshold between early-type and late-type hosts, 
as the latter (D/T > 0.575) demonstrate median populations of 
MATLAS-like dwarfs identical to the former. 

Dwarf galaxy populations in the MATLAS fields can also 
be described using luminosity functions (Fig. 3, central panel). 
Due to the large variance in the number of dwarfs per field and 
the discrete nature of dwarf counts, we characterised the isolated 
MATLAS luminosity function (red line) using the cumulative 
mean dwarf population brighter than a given M,. We note that 
2.8% of MATLAS dwarfs lack m, estimates from SEXTRACTOR 
or GALFIT, so we indicate the true mean dwarf population per 
isolated field with a red dotted line. As discussed in Sect. 3.1, 
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Fig. 3. Population of MATLAS-like dwarfs in TNG’s mock-observed fields is significantly lower than observed in the MATLAS low density fields. 
Left: dwarf counts in the 63’ by 69’ MATLAS fields (red) and MATLAS-like fields mock-observed in TNGSO (black). Dwarf galaxies that appear 
to be physically associated with the targeted hosts (by a radial velocity difference of <500 km s™!) are also shown in blue, while dwarfs bound 
to their host galaxy’s halo are drawn in green. Solid lines indicate fields around isolated MATLAS and MATLAS-like host galaxies that lack a 
companion with a magnitude of <Mx + 1 within 1 Mpc, while dotted lines represent the remaining fields. Dwarf populations of isolated TNG 
fields are inconsistent with these isolated MATLAS fields — even when including unassociated interlopers in the former — at a 70 confidence level. 
Centre: cumulative luminosity functions for MATLAS and TNG dwarfs in isolated fields. The solid line represents dwarf absolute magnitudes 
derived by assuming they lie at their supposed host galaxy’s distance (like the MATLAS dwarfs), while the dashed line adopts their true distances 
instead. Since around 3% of the MATLAS dwarfs lack m, estimates, the red dotted line indicates the true mean population of MATLAS dwarfs 
per field. Right: same as the left panel, but only counting dwarfs that lie within a square field with physical dimensions of 175 kpc centred upon 
each targeted host galaxy corresponding to the MATLAS field size at the near bound of 10 Mpc. 


M; is calculated for MATLAS dwarfs by assuming they lie at 
their presumed host galaxy’s distance, Dpost. Solid black and 
blue lines correspond to TNG dwarfs under the same assump- 
tion, while dashed lines show their true M, distribution. Solid 
and dashed lines are roughly equivalent for bound TNG dwarfs 
(green line) due to their physical proximity to their field’s cen- 
tral ETG and, thus, only one line is shown. To fully match TNG’s 
mean value of 0.74 bound satellites per isolated field, 93% of the 
MATLAS dwarfs would need to be interlopers unassociated with 
their supposed host. 

The one-degree fields in MATLAS cover a region around 
their central host galaxy that scales linearly with observation dis- 
tance — with a length ranging from 175 kpc at Dhost = 10 Mpc 
to nearly 800 kpc at the maximum host distance of 45 Mpc. To 
compensate for this dependence, we follow Habas et al. (2020)’s 
method and define a square window with dimensions of 175 kpc, 
corresponding to the minimum physical field size around each 
MATLAS ETG. The distribution of dwarfs located within this 
175 kpc window is plotted in the right-hand panel of Fig. 3. In 
contrast to the isolated MATLAS median of two dwarfs per win- 
dow, over 50% of isolated MATLAS-like fields in TNGSO con- 
tain no dwarfs within their 175 kpc window. According to the 
2-sample KS test, this discrepancy constitutes a 7.30 tension 
between the two distributions (D = 0.398, p = 3.1 x 10713), a 
significance even exceeding the tension over the full one-degree 
fields. 


3.3. Satellite galaxies and interloper contamination 


When comparing populations of dwarf galaxies in the MATLAS 
fields and TNG’s mock-observed analogues, the fraction of inter- 
loping background galaxies cannot be assumed to be similar a 
priori. The multi-stage visual inspection process in the MATLAS 
survey (Habas et al. 2020) can rule out a majority of highly lumi- 
nous non-dwarfs due to their distinct morphology and the pres- 
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ence of substructure. In TNG, however, we are limited to the use 
of MATLAS-like photometric cuts (Eqs. (3)-(6)) to attempt to 
achieve the same result. 

Habas et al. (2020) reported that 77-82% of MATLAS 
dwarfs with H1 velocities or distance estimates demonstrate 
radial velocities within 500kms~' of their fields’ central host 
galaxies. On the other hand, Heesters et al. (2023) recently 
obtained additional spectroscopy for 56 MATLAS dwarfs from 
VLT’s Multi Unit Spectroscopic Explorer (MUSE; Bacon et al. 
2010), finding that a much lower estimate of 57% is likely asso- 
ciated with their respective fields’ central ETGs using |A,| < 
500kms~!. In TNG’s mock-observed fields, a similar satellite 
fraction of 55% of all MATLAS-like dwarfs (43% in isolated 
fields) satisfies the same |A,| criterion, although the true fraction 
of dwarfs bound to their field’s central ETG is much lower at 
31% (16%). This counter-intuitive result of isolated fields pos- 
sessing lower association fractions is due to their targeted hosts 
tending towards lower halo masses (and hence smaller satellite 
populations); meanwhile, the number of background interlopers 
does not vary as strongly with the target host’s M200. 

We demonstrate further signs of a significant fraction of 
interlopers among the full sample of MATLAS dwarfs in Fig. 4, 
which plots the distribution of their distance-dependent proper- 
ties in gray. The MATLAS dwarfs are assumed to lie at their 
presumed host galaxy’s distance when calculating M, and Re. 
In TNG, the bound and background sample of MATLAS-like 
dwarfs follow similar intrinsic distributions in both M, and Re. 
However, the background dwarfs are fully distinct from the 
bound satellites when assuming they lie at a distance of Dhost. 
Due to the different distance ranges they inhabit, the intrinsic 
similarity between bound and background dwarfs also results in 
TNG’s two populations following distinct distributions in mg. W 
note that MATLAS dwarfs with distance estimates (black dotted 
lines) exhibit a significant excess at higher estimated luminosi- 
ties simply because brighter dwarfs are more likely to satisfy 
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Fig. 4. Properties of the MATLAS dwarfs and their simulated TNG analogues. TNG dwarfs bound to their host galaxy are plotted in black, while 
those unassociated with their supposed host are drawn in blue. MATLAS dwarfs are shown in shaded grey. Absolute magnitude, M,, and effective 
radius, Re, require distance estimates to derive, while apparent magnitude, mg, is a purely observational property. For the former two, estimates 
derived by assuming TNG dwarfs lie at their designated host galaxy’s distance are plotted in red. Of the MATLAS dwarfs, 15% have individual 
distance estimates. The M, and R, values calculated using these known distances are drawn as black dotted lines. The distribution of MATLAS 
dwarfs in M}, mg, and R, generally lies between the bound and background population in TNG. We also estimate the fraction of MATLAS dwarfs 
physically associated to their field’s host galaxy, fassoc, by averaging over TNG’s bound and background CDFs (green dashed line) with a ratio that 


minimizes the corresponding Kolmogorov-Smirnov statistic. 


the S/N > 5 requirement for their spectra (Habas et al. 2020) 
required to obtain sufficiently robust spectroscopic redshifts. 

Since the MATLAS dwarf catalogue should consist of both 
bound satellites and interlopers, we may consider its distribution 
in each of these three properties as a superposition of TNG’s 
mock-observed bound and background dwarfs at some given 
ratio. For each property, we identify the fraction of satellites 
associated with their presumed host, fassoc, Which minimises 
the summed CDF’s Kolmogorov-Smirnov statistic (green dashed 
line) with respect to the observed MATLAS distribution (black 
dashed line outlining the grey region). We obtained fassoc of 
0.55, 0.62, and 0.35 for M}, mg, and Re respectively. The opti- 
mised Re CDF does not trace the MATLAS distribution well, 
and its corresponding fassoe may be unreliable. Hence, we esti- 
mated that around 55—62% of all MATLAS dwarfs are associ- 
ated with their field’s central ETG. This fraction is consistent 
with Heesters et al. (2023)’s result of 57% (although we cau- 
tion that the association fraction of the brighter MATLAS dwarfs 
with available distance estimates may not directly translate to the 
full MATLAS sample). Conversely, our estimate is based on the 
full MATLAS dwarf sample, and remains valid if the photomet- 
ric properties of the MATLAS dwarfs are consistent with their 
TNG analogues. 


While we could attempt to alleviate the discrepancy in 
dwarf population between MATLAS and the mock-observed 
MATLAS-like fields in TNG by introducing a dominant frac- 
tion of unassociated background galaxies in the former, it seems 
likely that at least one-half of the observed MATLAS dwarfs 
are satellite galaxies of their field’s central ETG. Overall, TNG’s 
MATLAS-like fields demonstrate a comparable or greater degree 
of background contamination than MATLAS due to the lack of a 
visual inspection stage to remove massive interlopers with clear 
substructure. The tension between MATLAS and ACDM expec- 
tations in TNG50 remains robust at the >6c level. 


3.4. The impact of simulated resolution 


The dearth of MATLAS-like dwarfs in TNG50 (reported in 
Sect. 3.2) may be a direct result of an insufficient simulated 
resolution, as well as unreliable M} and Re estimates derived 
from too few stellar particles. In Fig. 5, we demonstrate that the 
baryonic resolution in TNG50 is sufficient to resolve a majority 
of MATLAS-like dwarfs and yield realistic intrinsic properties. 
Habas et al. (2020) estimates the distribution of stellar masses 
for 658 (30%) of MATLAS dwarfs with (g — i) colour estimates 
and clean photometry using Taylor et al. (2011)’s colour—mass 
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Fig. 5. Left: stellar masses of TNG subhalos by their association to their supposed host galaxy (black and blue lines) and estimates for the MATLAS 
dwarfs from (red line; Habas et al. 2020). The latter estimates were made by assuming that the observed dwarfs lie at the same distance as their 
field’s target host galaxy. Stellar masses for the subset of MATLAS dwarfs for which independent distance estimates are available are plotted as 
a red dotted line. TNG50’s baryonic resolution mgs is indicated by the hatched region. 97% of MATLAS-like dwarfs in TNG50 are sufficiently 
resolved with at least 100 stellar particles. Right: impact of simulation resolution on dwarf populations. The lower panel is taken from Fig. 2 and 
shows the distribution of luminous subhalo properties in R.-M, space in TNGSO (Mgas = 8.5 x 10* Mo), while the upper panel instead displays 
results for the lower-resolution TNG100 run (Mgas = 1.4 x 10° Mo). MATLAS-like dwarfs satisfying the selection criteria are plotted in black. The 
broadening in R, due to insufficient stellar particles occurs at around M, > —12 in TNG100, but only affects M, > —9 in TNG50 at a regime much 


fainter than the MATLAS-like dwarfs. 


relation (left panel, red line). The dotted red line shows the M, 
estimates re-computed for MATLAS dwarfs with independent 
distance estimates. The former consists of a fraction of back- 
ground galaxies that would systematically skew the distribu- 
tion towards lower M., and distances are only known for the 
most massive, easily observed dwarfs — thus skewing the latter 
distribution towards higher M,. As expected, we find the stel- 
lar masses of TNG’s mock-observed dwarfs between the two 
MATLAS CDFs. Finally, 97% of TNG50’s bound MATLAS- 
like dwarfs have a stellar mass above the equivalent of 100 stellar 
particles. 

In the right-hand panel, we show the intrinsic scaling rela- 
tion of all luminous subhalos in TNG50, as well as the lower- 
resolution TNG100 run (mga, = 1.4 X 10° Mo, mpm = 7.5 X 
10° Mo). Dwarfs considered to be MATLAS-like in each run 
are plotted as black points. In both runs, we identify a broad- 
ening in the Re distribution at higher M} (and lower M.) due 
to the difficulty in estimating the half-mass radius of a subhalo 
containing just a few star particles. In TNG100, this broadening 
begins at M, ~ 12 at stellar masses within an order of mag- 
nitude of most MATLAS-like dwarfs. This is not the case for 
TNGSO, wherein Re appears stable beyond M} < —10. We also 
recover a slight vertical shift in Re over the full subhalo distri- 
bution between runs, but this is likely an artefact of the larger 
force softening length adopted in TNG100 (185 pc) compared 
to TNG50 (74 pc). Hence, we argue that TNG50 has a baryonic 
resolution sufficient to resolve MATLAS dwarf analogues, and 
the discrepancy in its simulated population with the MATLAS 
fields is not a direct result of unreliable M, and (u4) estimates. 


4. Summary and conclusion 


The MATLAS low-to-moderate density fields survey low- 
surface brightness dwarf galaxies around massive early-type 
hosts (ETGs) at distances past 10 Mpc. Thus, they serve as an 
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excellent testbed for satellite galaxy abundances beyond our 
cosmic neighbourhood. To check whether dwarf populations in 
MATLAS are consistent with expectations from ACDM cos- 
mology, we mock-observed MATLAS-like fields in the high- 
resolution hydrodynamic simulation Illustris TNG-50 using pho- 
tometric selections faithful to the original survey. The simulated 
populations of MATLAS-like dwarfs sampled in this manner fol- 
low the standard scaling relations and demonstrate photometric 
properties within a range consistent with the MATLAS dwarfs, 
especially when assuming a 30—50% fraction of background 
interlopers unassociated with the targeted hosts. 

Strikingly, the simulated fields demonstrate significantly 
lower abundances of MATLAS-like dwarfs than found in the 
survey itself. Within isolated fields which only contain a single 
target ETG, MATLAS has a median of ten dwarfs with surface 
brightnesses of (ug) < 26.7 magarcsec~*, while TNG50 only 
has three such dwarfs. This result constitutes the basis for the 
“too-many-satellites” problem in observations at the >60 con- 
fidence level. We also account for the effect of host distance 
on the physical dimensions of the one-degree fields observed. 
We report that dwarf abundances within a central 175 kpc region 
around targeted ETGs demonstrate an similar tension at the 70 
confidence level. Furthermore, when using an viewing depth of 
100 Mpc, over 50% of mock-observed TNGSO dwarfs are back- 
ground interlopers that are not associated with their presumed 
host galaxy. On average, to match the simulated abundances of 
MATLAS-like dwarfs, 93% of the observed MATLAS dwarfs 
would need to be interlopers. This fraction is fully inconsistent 
with multiple estimates of the fraction of MATLAS dwarfs asso- 
ciated with their targeted ETG from the literature and this work, 
which do not drop below 50%. 

In this work, we consider whether this discrepancy in dwarf 
populations could have arisen from systematic biases. We find 
that the tension is unlikely to be a straightforward artefact of 
simulation resolution, since TNG50 sufficiently resolves 97% of 
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MATLAS-like dwarfs bound to their field’s targeted host with 
at least 100 stellar particles. The simulated dwarfs also follow 
the same scaling relations as found in MATLAS (Habas et al. 
2020), ruling out any significant systematic shift in their esti- 
mated effective radii. The effects of extinction from dust has 
not been considered, but this would further reduce the abun- 
dance of MATLAS-like dwarfs recovered in TNG50. We also 
do not account for the obfuscation of satellites in close projec- 
tion to the central ETG, but this would similarly work to reduce 
TNG50’s dwarf population and compound the problem. Finally, 
any incompleteness in the MATLAS fields would only exacer- 
bate the discrepancy by raising the true number of dwarfs, and 
the reported tension should be considered to be only a lower 
bound. 

Our findings lie in stark contrast with Carlsten et al. (2021), 
who reported the satellite luminosity functions of 12 Local 
Volume hosts to be fully consistent with expectations from 
dark matter-only cosmological simulations using a stellar-to- 
halo mass relation from Garrison-Kimmel et al. (2017) — despite 
covering a similar range in host mass and satellite surface bright- 
ness. We point out, however, that we did not assess the full 
luminosity function of satellites around the ATLAS*> ETGs, 
but instead we imposed strict photometric selection cuts to 
match our mock-observed fields with the MATLAS observa- 
tions. This resulted in a comparison between better defined 
sub-populations of dwarf galaxies. Furthermore, Carlsten et al. 
(2021) adopted the lower-resolution TNG100 run as their fidu- 
cial hydrodynamic simulation, which struggles to resolve dwarfs 
with luminosities corresponding to the fainter MATLAS dwarfs. 
It is also interesting to note that the range of stellar masses 
in which Miiller et al. (2024) reported the luminosity func- 
tion of M 83 to be most discrepant with theoretical predic- 
tions (107—108 Mo) coincides with the majority of MATLAS- 
like TNG50 dwarfs in this work. Our results suggest that mod- 
els of galaxy formation within the ACDM framework may 
still struggle to reliably reproduce dwarf abundances beyond 
the Local Volume. Thus, the too-many-satellites problem is 
extended from M 83 to across 150 fields in the MATLAS 
catalogue. 
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Appendix A: Impact of host isolation 


We verified whether the results obtained in this work is sensi- 
tive to our adopted isolation criterion of 1 Mpc imposed on the 
host galaxies mock-observed in TNG50. In Fig. A.1, we plot the 
mean dwarf counts in MATLAS-like fields as a function of their 
hosts’ distance to their closest companion. Until around 600 kpc, 
where the satellites of the companion host begin to encroach 
on the target’s dwarf population, the number of MATLAS-like 
dwarfs within the mock-observed fields (especially those phys- 
ically associated to their presumed host galaxy) remains rela- 
tively stable. We conclude that our specific choice of 1 Mpc in 
our isolation criterion is unlikely to significantly contribute to the 
observed discrepancy in dwarf populations. Imposing a stricter 
degree of isolation would only serve to artificially reduce the sig- 
nificance of the discrepancy due to the smaller sample of MAT- 
LAS fields. 
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Fig. A.1. Mean dwarf populations of isolated MATLAS-like TNG50 
fields, binned by the 3D distance between target hosts and their closest 
bright companion with a magnitude of Mx + 1. Black, blue, and green 
points correspond to samples of all MATLAS-like dwarfs, those with 
a radial velocity within 500 kms! of their presumed host, and dwarfs 
considered to be gravitationally bound to their host’s halo in Subfind — 
error bars in their respective colours denote their 1o spread. The upper 
panel shows the cumulative distribution of mock-observed MATLAS- 
like field realisations. This work’s adopted isolation criterion of 1 Mpc 
is indicated by the grey dashed line. 


